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ABSTRACT 

This  paper  considers  the  flow  past  a  flat  plate.  When  the  flow  is 
uniform  the  Blasius  boundary  layer  forms  on  the  plate.  Herein  the  usual 
uniform  free  stream  is  perturbed  by  a  small  amplitude  field  of  random  three- 
dimensional  turbulence.  Using  a  numerical  technique  the  statistical  response 
in  the  boundary  layer  due  to  the  turbulence  in  the  free  stream  is  computed  in 
terms  of  auto-  and  cross-correlation  functions.  The  field  variables  are 
expanded  in  finite  Chebyshev  series  and  the  space-time  correlations  of  these 
functions  are  found  by  forming  correlation  matrices  with  the  vectors  of  random 
coefficients  of  the  Chebyshev  expansions. 

The  main  result  of  the  evaluations  is  the  fact  that  three  dimensional 
random  infinitesimal  free  stream  turbulence  produces  a  fully  three  dimensional 
response  in  the  boundary  layer.  The  excited  field  includes  the  presence  of 
streamwise  vorticity  as  well  which  was  previously  thought  to  be  due  to 
secondary  instability.  /T\ 
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A  STATISTICAL  ANALYSIS  OF  THE  EFFECT  OF  FREESTREAM 
TURBULENCE  OF  THE  BLASIUS  BOUNDARY  LAYER 
Thomas  J.  Bridges  and  Philip  J.  Morris 

1.  Introduction 

The  subject  of  this  paper  is  an  analysis  of  the  response  of  a  boundary  layer  (the 
Blasius  boundary  layer)  where  the  free  stream  contains  random  three-dimensional  turbulence 
of  small  amplitude.  The  statistics  of  the  free  stream  turbulence  are  considered  known, 
and  from  these  and  the  governing  equations  a  statistical  solution  is  found  in  the  boundary 
layer. 

Previus  theoretical  work  on  this  problem  up  until  1975  has  been  summarized  by  Rogler 
and  Reshotko  (1975)  and  this  may  be  found  in  their  paper.  In  their  work  Rogler  and 
Reshotko  (1974,  1975)  studied  the  problem  in  question  when  the  free  stream  disturbances 
are  deterministic  and  two-dimensional.  The  free  stream  is  composed  of  a  uniform  flow  plus 
a  set  of  inviscid  point  vortices  which  convect  with  the  mean  flow.  The  analysis  is 
limited  by  the  two-dimensionality  but  their  results  show  a  strong  interaction  between  the 
boundary  layer  and  the  free  stream  with  excited  amplitudes  of  the  boundary  layer  vorticity 
reaching  twice  the  free  stream  value.  Their  result  also  shows  that  there  is  a  correlation 
between  the  length  scale  of  the  point  vortices  and  the  distance  from  the  leading  edge. 

In  this  paper  the  work  of  Rogler  and  Reshotko  is  extended  to  include  a  three 
dimensional  analysis  and  the  free  stream  disturbances  are  taken  to  be  a  field  of  random 
three  dimensional  turbulence.  This  free  stream  turbulence  is  three-dimensional  in  space, 
is  time  dependent,  and  is  composed  of  a  fluctuating  velocity  and  pressure  field  of  small 
amplitude.  The  properties  and  characteristics  of  the  random  free  stream  turbulence  are 
identified  by  specifying  a  multi-dimensional  autocorrelation  function.  This 
autocorrelation  function  includes  length  and  time  scales  as  well  as  spectral  information. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 .  This  material  is 
based  upon  work  supported  by  the  National  Science  Foundation  under  Grant  No.  DMS-8210950, 
Mod.  1. 
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The  response  is  composed  of  the  Blasius  solution  plus  a  perturbation  in  the  boundary 
layer  due  to  the  presence  of  external  free  stream  turbulence.  Since  the  forcing  function 
is  random,  the  excited  response  will  also  be  random.  Therefore,  the  solution  is 
determined  in  terms  of  auto-  and  cross-correlations  of  the  excited  velocity,  pressure,  and 
vorticity  fields. 

This  analysis  can  be  considered  a  natural  extension  of  the  theory  of  hydrodynamic 
stability.  In  hydrodynamic  stability  theory,  qualitative  solutions  are  obtained  for  a 
wide  range  of  parameters.  By  treating  the  boundary  layer  as  a  multi-degree-of-f reedom 
vibrating  system,  the  proper  linear  combination  of  these  natural  modes  of  the  boundary 
layer  can  be  summed.  This  results  in  a  statistical  solution  in  the  boundary  layer  which 
is  due  to  a  particular  random  forcing  function  in  the  free  stream. 

Before  proceeding  to  the  details  of  the  analysis  of  the  three-dimensional  boundary 
layer,  the  background  for  the  analysis  can  be  illustrated  by  considering  the  simple 
analogous  problem  of  a  linear  SDOF  system  under  stochastic  loading.  Such  a  system  has  the 
following  governing  equation  and  boundary  condi tons.. 

L<p  =  q(t)  (1.1) 

„{0)  «  <pM0)  -  0  (1.2) 


=  +  28u„ 


8  is  a  damping  parameter,  and  idq  is  the  undamped  natural  frequency. 

In  the  first  analysis,  the  homogeneous  operator  is  studied.  The  eigenproperties  can 
be  used  to  form  a  transfer  function  which  is  used  utimately  to  construct  the  complete 
solution.  By  taking  the  Fourier  transform  of  (1.3),  the  result  is 
2  2  * 
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The  statistical  properties  of  the  random  forcing  function  must  be  defined  based  on 


some  physical  observation.  Taking  q(t)  to  be  statistically  stationary  an 


autocorrelation  can  be  defined. 


R  (t)  =  <q*(t),q(t  +  t)> 

qq 


where  the  brackets  denote  an  ensemble  average  and  *  denotes  complex  conjugate. 


Alternatively  a  spectral  density  function  S^(uj)  may  be  defined,  which  is  related  to 


R  (t)  by  the  Fourier  transform.  If 

qq 


q(t)  *  /  q(m)e  iutdui 


<q*  ( €01 ) ,  q(  0J2  )  >  =  S  ( t»>1 )  6  ( o>2  -  u^) 


where  the  delta  function  arises  because  q(t)  is  stationary. 


The  next  step  is  to  obtain  a  solution  for  e(t)  in  time.  Taking  a  Fourier  transform 


of  Lp  =  q  and  defining 


F(w)  =  [— u>2  -  2i8u0<ii  +  “f.2  ]  1 


then  the  solution  is 


n>(t)  ■  /  F(ui)q(to)e  utda) 


(1.10) 


However  only  the  statistical  properties  of  q(w)  are  known.  Therefore  an 


autocorrelation  function  for  <p ( t )  must  be  obtained. 


R  ( x )  =  <<p*(t)  ,9(t  +  t)> 
<P<P 


(1.11) 


Substituting  (1.10)  into  (1.11)  and  using  (1.8)  the  result  is 


R  (t)  =  /  |  F ( u> >  |  S  (u)e  110  dw 

<P9  qq 


(1.12) 


Once  S  (m)  is  defined  explicitly  the  integral  can  be  inverted  using  residue  theory, 

qq 


where  the  poles  are  those  given  in  (1.5). 


Although  the  addition  of  three  spatial  dimensions  in  the  boundary  layer  problem  will 


add  considerable  complexity  to  the  analysis,  the  basic  approach  is  the  five  steps  given 
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above:  (a)  the  equations  governing  the  boundary  layer  are  derived)  (b)  the  eigenvalues 
and  eigenfunctions  of  the  operators  are  analyzed  which  are  used  to  construct  a  transfer 
function;  (c)  the  statistical  properties  of  the  random  free  stream  turbulence  are  defined 
by  recourse  to  experimental  observations;  (d)  the  solution  is  found  in  real  space;  and  (e) 
the  multi-dimensional  statistics  of  the  excited  velocity  and  vorticity  field  in  the 
boundary  layer  are  obtained.  These  results  can  be  further  distilled  to  determine  such 
parameters  as  excited  length  scales,  and  the  transverse  distribution  of  the  excited 
response  in  the  boundary  layer. 

In  Section  2  the  equations  governing  the  flow  past  the  flat  plate  are  derived.  The 
analysis  of  the  eigenvalues  and  eigenfunctions  of  the  derived  operators  have  been  reported 
in  a  companion  paper  (Bridges  and  Morris,  1985).  In  Section  3  the  free  stream  turbulence 
is  analyzed  and  the  necessary  multidimensional  autocorrelation  function  for  the  free 
stream  turbulence  is  specified.  In  Section  4  the  inhomogeneous  equations  are  solved  for 
the  velocities,  pressure,  and  vorticities.  Fourier  transforms  are  used  in  the 
x,  z,  and  t  dimensions  for  obtaining  the  solution,  and  the  solution  in  the  y-dimension 
is  expressed  in  terms  of  a  series  of  Chebyshev  polynomials.  The  y  dimension  which  goes 
from  0  to  «  is  transformed,  for  the  Chebyshev  expansion,  to  Ye  t-1,+1]  using  the 
transformation 

Y  »  (y  -  2)/<y  +  2)  (1.13) 

with  metric 

m(Y)  -  dY/dy  =  <1  -  Y)2/4  (1.14) 

In  Section  5  the  statistics  of  the  excited  field  are  derived.  The  multi-dimensional 
autocorrelation  functions  are  expressed  as  double  Chebyshev  expansions,  and  convenient 
expressions  for  the  associated  matrix  of  coefficients  are  constructed. 

Residue  theory  and  numerical  integration  are  used  to  invert  the  Fourier  integrals 
which  result  in  the  correlation  functions.  The  details  necessary  to  complete  this  are 
described  in  Section  6.  In  Section  7  the  correlations  are  computed  for  R  =  1000,  2000, 

-4— 


and  3000  and  results  are  shown  of  the  distribution  across  the  boundary  layer  of  the  mean 
square  velocities  and  vorticities  and  the  Reynolds  stress .  By  evaluating  the  expressions 
for  various  values  of  the  spanwise  separation  typical  spanwise  length  scales  of  the 
excited  solution  are  determined.  The  main  result  is  the  prevalence  of  three- 
dimensionality  in  the  solutions.  It  is  shown  that  if  the  random  free  stream  turbulence  is 
three-dimensional,  then  the  excited  solution  will  be  patently  three-dimensional. 

Comparison  with  the  experimental  results  of  Klebanoff,  et  al.  (1962)  shows  excellent 
agreement. 


*>: 


-5- 


Accosslon  Vor 

NTTS  GRA&l 
DTIC  TAB 
Ur.cruiowiced 
Ju:;l  li'  loot  ion_ 


□ 

□ 


By - 

JJiatribution/ 

Availability  Codes 
Avail  and/or 


Dlst 


Special 


(2.1) 


The  equations  governing  the  flow  field  are  the  Navier-Stokes  equations 
V  •  U  «  0 

—  +  (u  •  v)u  +  -  Vp  -  vAU  =  0  (2.2) 

ot  p 

where 

U  -  Ui  +  Vj  +  Wk 

A  is  the  Laplacian,  and  the  boundary  conditions  are 

0  «  V  »  W  »  0,  at  y  =  0,  for  all  z  and  for  x  >  0 


0  ♦  u  +  eu 
«  s 

V  *  eu 

s 

W  >  ew 

s 

P  *  p  +  ep 
“  *s 


as  y  +  “ 


(2.3) 


Since  the  free  stream  turbulence  has  amplitude  e,  a  small  parameter,  it  is  proposed  that 

the  response  in  the  boundary  layer  can  be  expressed  as  a  power  series  in  e ; 

2 

U(x,y,z,t)  -  uo<x,y,z)  +  eu^x.yjZ.t)  +  0(e  ) 

2 

V(x,y,z,t)  »  vQ(x,y,z)  +  ev^x.y.Zjt)  +0(e  ) 

2 

W(x,y,z,t)  -  wQ(x,y,z)  +  ew^x.y.z.t)  +  0(e  ) 

2 

P(x,y,z,t)  -  PQ(x,y,z)  +  ep1(x,y,z,t)  +  0(e  )  (2.4) 

The  expansions  (2.4)  may  be  substituted  into  (2.1),  (2.2),  and  (2.3)  and  terms 
proportional  to  like  powers  of  c,  are  equated  to  zero.  Only  the  zeroth  and  first  order 
problems  are  considered.  The  zeroth  order  problem  reduces  to  the  Blasius  equation  if  the 
usual  assumptions  of  boundary  layer  theory  are  used.  Then  retaining  the  parallel  flow 
approximation,  usual  in  hydrodynamic  stability,  the  first  order  problem,  after  rendering 
the  variables  dimensionless  using  a  boundary  layer  thickness  and  the  free  stream  velocity, 
and  shifting  the  dependent  variables  to  render  the  far  field  boundary  conditions 


homogeneous,  i6 


A P'  =  -2  —  {v*  +  v  ) 

3y  dx  s 


,  32u  3v 

L  v*  =  (1  -  u  )  Av  +  - ^  — - 

4  0  3x  s  .  2  3x 

3y 


V  =  - 


3p'  3us  ^0 

>'  -  +  (1  -  u„)  — r—  -  v  - 

3x  0  3x  s  3y 


lr+  M 


where  the  L2  and  L4  operators  are  defined  as 


t  _  r  9  ,  5  ^  a  1 

V  =  Ljt  +  U0  3x  ~  R 


V  =  L2i» 


(2.10) 


and  the  boundary  conditions  are 


u'(x,0,z,t)  «  -ug(x,0,z,t) 


v' ( x, 0,z, t)  =  -v  (x,0,z,t) 


w'(x,0,z,t)  =  -w  <x,0,z,t) 


p' ( x , 0 , z , t )  =  -ps(x,0,z,t) 


u'.v'.w'jp'  +0  as  y  + 


(2.11) 


where  the  primed  variables  are  differences  between  the  first  order  variables  and  the  far- 


field  boundary  values,  e.g.  v’  =  v1  -  vg,  etc. 


The  operator  is  the  well  known  Orr-Sommerf eld  operator.  The  eigenvalues  and 


eigenfunctions  of  this  operator  play  an  important  part  in  the  construction  of  solutions  to 


the  randomly  forced  problem  at  hand.  A  study  of  the  necessary  aspects  of  this  operator, 


as  well  as  the  Lj  operator,  have  been  detailed  previously  by  Bridges  and  Morris  (1985). 


In  the  next  section  the  equations  for  the  free  stream  turbulence  are  analyzed,  and 


their  statistical  properties  are  formulated. 
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In  this  section,  aspects  of  the  forcing  function,  the  free  stream  turbulence,  are 


defined.  Using  the  governing  equations  for  the  free  stream  turbulence,  a  relationship 
between  the  transverse  velocity  and  the  other  velocities  is  derived.  Therefore  only  the 
statistical  properties  of  the  transverse  velocity  need  be  specified.  An  autocorrelation 
function  for  the  transverse  velocity  is  then  specified  based  upon  known  experimental  data 
gathered  in  wind  tunnels  on  decaying  turbulence.  The  other  statistical  properties  of  the 
free  stream  turbulence  can  be  derived  from  this  function. 

The  equations  which  govern  the  free  stream  turbulence  are  the  first  order  equations, 
derived  in  the  last  section,  evaluated  in  the  free  stream  as  y  +  “.  The  equations  can  be 
solved  for  a  relationship  between  ua,  vg,  and  wg.  Defining  transform  mates  by 


v  (at y;8,w) 
s 


1  rtr  ,  -i( ax+Sz-ut)  .. 

- -  Iff  v  (x,y,z,t)e  dxdzdt 

( 2ir)  J  —  8 


vg(x.y.z.t) 


///  v  (ajy ; 8 ,u)e 


i(ax+Bz-iot) 


dadgdu 


then  the  other  velocity  components  are, 


u  <x,y,z,t)  -  fff  — — —  vg(a,y,8,u>)e 
y  -«  a  +8 


dadgdw 


w  ( x,y,z,t)  -  -jp  fff  18  v  (a,y,8,io)ei(aX+®Z  “^dadBdm 
3  y  -»  0^+8 


(3.3) 

(3.4) 


which  may  now  be  used  to  obtain  statistics  when  the  correlation  function  for  vg  is 
specified. 

An  autocorrelation  function  with  sufficient  generality  is  chosen  to  represent  the 
statistical  properties  of  the  transverse  velocity.  As  most  of  the  experimental  data  on 
the  properties  of  small  amplitude  free  stream  turbulence  is  taken  in  wind  tunnels,  the 
qualitative  features  of  this  flow  are  used  as  a  guide.  The  autocorrelation  function  is 


chosen  to  be  stationary  in  time,  homogeneous  in  the  transverse  and  spanwise  directions, 
and  inhomogeneous  (decaying  turbulence)  in  the  streamwise  direction.  To  represent  the 
homogeneous  aspects  of  the  correlation  function  a  general  Gaussian  form  is  chosen,  which 
includes  a  specified  integral  length  or  time  scale.  For  example,  in  the  transverse 
direction  the  correlation  function  is  taken  to  be 


Ry(y1/y2)  *  e*p<-”<y2  -  y//4ij>  o.S) 

where  l  is  the  integral  length  of  the  turbulence  in  the  transverse  direction  defined  by 

00 

i  •  (  expf--n(~- )  )dn  (3.5a) 

y  0  y 

where  n  =  y2  "  y.)-  This  length  scale  is  representative  of  the  characteristics  of  the 
largest  eddies  in  the  turbule--  In  the  streamwise  direction  a  decaying  Gaussian 
function  is  used, 


R 

x 


(x1,x2) 


A(x  ^ )  exp(-s(x2 


(3.6) 


In  wind  tunnel  experiments  on  decaying  turbulence,  while  the  intensity  of  the  turbulence 
is  decaying,  the  length  scale  is  growing.  However,  most  experiments,  for  example 

those  of  Compte-Bellot  and  Corrsin  (1958),  suggest  that  tx  grows  like  the  square  root  of 
the  distance  downstream,  and  since  it  is  here  rendered  dimensionless  by  a  boundary  layer 
thickness,  which  also  grows  like  the  square  root  of  the  distance  downstream  and  is 
considered  to  be  a  local  constant,  in  congruence  with  the  parallel  flow  approximation, 

can  also  be  considered  as  a  local  constant.  A  plausible  expression  for  the  decay  of 
the  amplitude  of  the  correlation  would  be  an  exponential,  A(x.,)  =  AQexp(— yx^.  Later  in 
the  calculations,  the  initial  amplitude  can  be  chosen  to  be  unity  (essentially  included  in 
the  small  parameter  e),  and  the  results  scaled  accordingly. 

The  fact  that  the  turbulence  field  is  superimposed  on  a  uniform  velocity  field 
suggests  that  the  proper  correlation  function  should  be  a  convected  correlation.  Such  a 
convected  field  is  distinguished  by  the  tendency  for  maximum  correlation  to  lie  along  a 


space-time  trajectory  where  (x2  -  x^ )  -  Uc(t2  -  tj)  is  zero.  Consequently,  since  the 


turbulence  is  stationary,  and  the  convection  velocity  is  unity,  the  streamwise  correlation 
is 

2  2 

Rx*X1'X2,T*  *  A*xi)  exp(-ii(x2  -  x,-  t)  /4lx  )  •  (3.7) 

Therefore  the  complete  autocorrelation  function  for  the  free  stream  transverse  velocity  is 

WX2'y1'y2*C'x)  ' 

.2  2  2 


(x,  -  x  -  t)  ... 

Mx1)*xP  [-  - “1 -  +  ~~2~  +  + 

l  l  i  t 

x  y  z  t 


(3.8) 


where  £  =  z2  -  z^. 


Much  of  the  later  analysis  will  take  place  under  the  cloak  of  Fourier  integrals. 
Therefore  the  Fourier  transform  of  the  x,z,t  portion  of  the  autocorrelation  function  is 
necessary.  Splitting  Rs  into  two  parts, 

WX2'yl'y2'C'T)  "  Ry(y  l'y2)C(xrX2'C'T> 

and  defining  the  four  dimensional  Fourier  transform  as. 


(3.9) 


.  .  •  -i(a  x  +a  x  +8C“iot) 

Qfo^.o^.B.io)  -  - -  Iff!  QCx^x^jC.tle  dxldx2dCdT 

(2w)  -» 


(3.10) 


then  the  desired  result  is 


e<0‘ra2,8'“> 


At  l  t 
0  x  z  t 


2w 


r - ’ - 1 

4  ly  +  i(a1  +  a2 ) J 


exp[-0  /x] 


(3.11a) 


where 


0A 


l2*  +  t2B2  +  (u  -  a  )2l2 

x  2  z  2  t 


(3.  1 1b) 


In  the  subsequent  analysis  the  solution  in  the  transverse  direction  will  be  obtained 
by  finite  expansion  of  the  solution  in  a  Chebyshev  series.  This  requires  that  the  forcing 


function  also  be  expanded  in  a  Chebyshev  series, 

N 


vg ( x,y ,z, t)  -  Z'  qn(x,z,t)Tn<¥) 
n-0 


(3.12) 
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where  the  prime  signifies  that  the  first  term  is  to  be  halved,  and  y  has  been 
transformed  to  Y(y),  where  Y  e  [-1,+1].  It  should  be  pointed  out  here  that  no  attempt 
is  made  to  actually  expand  the  free  stream  turbulence  as  a  deterministic  function.  The 
qn(x,z,t)  are  random  amplitudes.  Consequently ,  only  statistical  properties  of  qn  will 
be  used.  Therefore,  the  accuracy  of  the  Chebyshev  representation  of  the  free  stream 
turbulence  will  rest  on  how  well  the  autocorrelation  of  the  free  stream  turbulence  can  be 
expanded  in  a  double  Chebyshev  series;  and  this  can  be  done  quite  accurately. 

In  terms  of  vg,  the  autocorrelation  function  Rg  is 

RS(x1'x2'y1'y2'C'T)  =  <v*(x,<y1'z't>vs(x2,y2,Z  +  ;,t  +  T»  (3.13) 

where  *  denotes  the  complex  conjugate  and  the  brackets  denote  the  ensemble  average. 
Substituting  (3.12)  into  (3.13)  results  in 

<vs(x1'y1'z't,Vs(X2'y2'Z  +  C't  +  T)>  = 

N  N 

y*  V1  <q^(Xl  .z.tlq^tXj  ,z  +  c,t  +  T )  >  TjjtY,  )Tm(Y2)  (3.14) 

n=0  m=0 


Due  to  the  orthogonality  of  the  Chebyshev  polynomials,  the  auto- 
the  qn  are 

<lCX1'Z't)<*m<X2'X  +  C>t  +  T)>  =  C(X1'X2'C'T)Rmn 


and  cross-correlations 


(3. 15) 


of 


where 


and 


where 


Q(x1 ,x2 ,;,t)  =  A( x^ )  exp 


.  1  T  (Y.)dY,  1  K  (Y,,Y,)T  (Y,)dY_ 

-  _5.  (■  n  1  1  ,  y  12m  2  2 

mn  0  '  -  - J 


/I  -  Y,2 


-1 


/i-yT2 


(3.16) 


(3.17) 


R  ( Y  ,  Y  )  =  exp  [  - 


[y,<V  -  y.(Y,)]2i  • 


(3.16) 


This  bears  out  the  previous  conjecture,  that  the  accurate  representation  of  the  random 


transverse  velocity  in  the  free  stream  depends  on  how  accurate  the  y-  autocorrelation 
function  can  be  expanded  in  a  double  Chebyshev  series. 

As  noted,  the  x,  z,  t  behavior  will  be  analyzed  under  Fourier  Integrals. 
Therefore  the  necessary  analogy  of  equation  (3.15)  in  the  Fourier  domain  is 


R  Q(-a.,a  ,  )6(ui  -  to, )  S  ( 6-  -  8,) 

run  izii  4  l  4  i 


(3.19) 


where  5(x)  is  the  Dirac  delta  function,  is  given  by  (3.17),  and  Q  is  given  by 

(3.11). 

This  completes  the  specification  of  the  statistical  properties  of  the  forcing 
function.  It  is  a  function  of  five  parameters;  the  integral  length  scales  in  the  x,  y, 
and  z  directions,  the  integral  time  scale,  and  the  decay  rate  of  the  turbulence  in  the 
streamwise  direction.  In  the  next  section  a  method  for  determining  expressions  for  the 
excited  field  in  the  boundary  layer  is  derived. 


Section  4.  Solution  for  the  First  Order  Field 

The  solution  for  the  first  order  excited  field  in  the  boundary  layer  follows  the 
analysis  of  the  homogeneous  problem  outlined  in  Bridges  and  Morris  (1985).  The  non¬ 
constant  coefficients  in  the  equations  are  function  of  y  only.  Therefore  a  Fourier 
transform  is  used  in  the  x,  z ,  and  t  dimensions.  The  y  dimension,  which  extends 
from  0  to  is  transformed  to  the  Chebyshev  domain,  V  e  1-1, +1],  using  the 

algebraic  transformation, 

Y  -  <y  -  2)/(y  «■  2)  (4.1) 

with  metric 

m<y)  -  dY/dy  -  (1  -  Y)2/4  (4.2) 

The  y  dependence  of  the  solutions  is  then  expressed  in  terms  of  a  finite  Chebyshev 
series. 

The  steps  necessary  to  obtain  an  expression  for  the  transverse  velocity  will  be  gone 
through  in  some  detail  for  illustrative  purposes  and  then  the  solutions  for  the  remainder 
of  the  field  variables  will  be  only  briefly  described.  Further  details  can  be  found  in 
Bridges  (1984). 

To  render  the  boundary  conditions  on  the  transverse  velocity  completely  homogeneous, 
the  function  v'  defined  in  (2.13)  is  split  into  two  parts, 

vi  (x,y,z,t)  =  vs(x,y,z,t)  -  g.((x,y,z,t)  +  v(x,y,z,t)  (4.3) 

where 

3v 

g^x.y.z.t)  =  va(x,0,z,t)e'y  +  (vs(x,0,z,t)  +  -^  (x,0,z,t)  )ye  y  (4.4) 

and  v,  the  only  unknown  function,  is  governed  by  an  inhomogeneous  partial  differential 
equation  with  homogeneous  boundary  conditions.  Taking  a  Fourier  transform  of  this 
equation  in  x,  z,  and  t,  where  the  Fourier  transform  mates  are  defined  in  equation 
(3.2),  and  multiplying  the  equation  by  -R  results  in  the  governing  equation  for  v, 

2  d2u 

L  v  =  iaR(u.  -  1)[— rr  -  (a2  +  S2)]v  -  iaR  — v  +  L  g  (4.5) 

4  0  .2  4  S  4S41 

dy  dy 
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£ 


where 


d2u„ 


2  2 

=  [iR(<o  -  ou0)  +  -(a2  +  e2)]f-^  -  (a2  +  e2)]i  +  iaR - j  ' 


(4.6) 


dy  dy  dy 

is  the  well  known  Orr-Sonmerfeld  operator.  The  equation  is  transformed  to  the  Y  domain 

A 

using  the  algebraic  transformation  (4.1),  and  the  unknown  function  v,  and  the  known 


functions  v  and  uQ(Y)  are  expanded  in  finite  Chebyshev  series. 


-  N  . 

v(a,y,B,w)  -  T'  v  (a , 8 ,(o)T  ( Y) 
n  n 


(4.7) 


n“0 


Using  the  Chebyshev  discretization  approach  described  in  Bridges  and  Horris  (1985),  the 


equation  is  integrated,  the  Chebyshev  expansions  for  v,  vg,  and  uQ  are  substituted 


into  the  equation,  along  with  the  boundary  conditions,  and  the  result  is  a  matrix 
equation. 


[D4(o,8,u)]{vn}  -  [H0(a,8,u)]{qn} 


(4.8) 


where  the  Reynolds  number  is  taken  to  be  specified,  and 


D  (a, 6,(0)  -  C  <B,w)a4  +  C,(8,(o)o3  +  C  (B,co)a2  + 
4  0  12 


«  (C(a, 6, util  (H-<a, 8, all  „ 

-  tit  l - SuX Si - -  *  11  -  =,>]<V* 


dadBdco 


(4.16) 


If  the  expression  for  {q^  were  Known,  the  integral  (4.16)  could  be  inverted  using 
residue  theory  where  the  poles  would  arise  from  A(a,B,k>)  =  0,  and  from  any  poles  of 
(q^.  In  a  subsequent  section  it  will  be  shown  that  the  poles  of  A(a,B»w)  are 
exponentially  larger  than  the  poles  due  to  the  forcing  function.  Therefore,  in  the 
calculations  the  contribution  of  the  matrix 
(I  -  Gj]  to  the  overall  solution  is  negligible. 

The  pressure  p'(x,y,z,t)  satisfies  the  Poisson  equation  given  in  (2.14).  Using  the 
y-momentum  equation,  it  can  be  shown  that 


-15- 


a  +  B 


2  2 
a  +  8 


lia<1  “  un>  “3T  +  ia  VJ 


(4.17) 


where  v'(x,y,z,t)  =  v(x,y,z,t)  -  g^x.y.z.t).  The  solution  for  p1  in  real  space  is 


p^x.y.z.t)  -  pg(x,y,z,t)  +  ///  p'(a,y,8,u)e 


i(ax+8z-u>t) 


dadBdui  • 


(4.18) 


The  governing  equation  for  the  streamwise  velocity  u* (x,y,z,t)  is  given  by  equation 
(2.16).  If  a  boundary  condition  function  is  used  to  render  the  boundary  conditions 


homogeneous,  then  u^  has  the  expression 

u^x.y.z.t)  =  u(x,y,z,t)  +  ug(x,y,z,t)  -  g2(x,y,z,t) 


(4.19) 


with  g2  ■  e-y  ug(x,0,z,t),  and  the  only  unknown  is  u(x,y,z,t).  After  substitution  of 
(4.19),  and  subsequent  exercise  of  the  Fourier  transform  in  x,z,  and  t,  the  governing 
equation  for  u  is 


V  ‘  “i<XR<1  "  U0,Us  +  R  dty(V'  +  V8>  *  iORp’  +  L292 


2*  a 
L2u  »  +  [iR(w  -  auQ)  -  (a2  +  B2)]u  . 

dy 

A  * 

However  u  is  related  to  v  by 
s  s 

dv 

io  8 


3  a2  +  B2 

and  p'  can  be  expressed  in  terms  of  v  and  vg.  Then  (4.20)  becomes 

*  *  *  *  B2R  dU0  *  82R  du0  *, 

V  -  L292  +  a2  +  62  dy  vs  +  a2  +  g2  dy  V 

+  -r2-!  ^’+  :iR(“  -  aV  -  <°2  + 

a  +  8  dy 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


vV/V.vs/v- 


**,  f  H  •  •'  *  4  **  ••  *’  « 

-  •  -j- 


\  -\  v\  ■,  \  ■  •  / 


the  unknown  u  is  then  expended  in  a  finite  Chebyshev  series 


u(aiy;6,u>)  -  V'  u  (a,8,w)  T  (Y) 
N-0  n  n 


(4.24) 


Substituting  (4.24)  as  well  as  the  Chebyshev  series  for  the  free  stream  velocity,  and  the 
known  transverse  velocity,  the  following  matrix  equation  results, 

AAA 

[D2(a,8,u)l (un)  -  (H1 (a,8,u>) (vn)  +  [H3(a,6,w>] {q^J  (4.25) 


where  {v  }  is  given  by  (4.8).  The  poles  of  0.(a,8,ci>)  are  all  quite  stable  (a  proof  of 


this  is  given  in  Appendix  C)  and  will  lend  little  to  the  overall  solution  by  comparison  to 


the  poles  of  {v^}.  Therefore, 


[D  (a, 6,0))]  [H  <a,8,<*>>]  (C(a,B,u)]  [H0(a,8,u)] 


<«n> 


(4.26) 


A(  o,B ,u) 

Taking  into  account  the  expressions  for  ug  and  g2 ,  the  first  order  streamwise  velocity 
has  a  Chebyshev  expansion. 


u. (x,y,* ,t)  -  y>  u.  (x,s,t)  T  (Y) 
N-0  "  n 


(4.27) 


with  the  vector  {u,nJ  given  by 


»  (D,  ]  [H  ]  [C]  [H  ] 
<*1n>  *  ///(  2"(o-78,-M-  — •  * 


-2-  . . 2  U  -  G2]  ID,]  ]{qn}ei<e,X+Sl'“t)dad8da)  (4.28) 

a  +  fT 

and  the  matrices  G2  and  D,  are  described  in  Appendix  A. 

Using  a  similar  procedure,  the  complete  expression  for  w,,  the  first  order  spanwise 
velocity  is 

N 

w,(x,y,z,t)  -  J*  w,n(x,z,t)  Tn(Y)  (4.29) 

with  the  associated  vector  given  by 
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V, 


i 


{w  }  -  /rr r _ l - L 

'"in*  1>n  A(a,8,w) 


i8(D  ']  [H  J  [C]  (H  ] 

-TO -  +  7T7  11  -G2)(D11](qn>  • 

a  +  p 


ei<ax+6z-“t)dad0du 


(4.30) 


The  excited  first  order  vorticity  field  in  the  boundary  layer  is  also  of  considerable 
interest.  The  streamwise  vorticity  is  given  here.  The  other  two  components  can  be 
derived  in  a  similar  manner.  The  first  order  streamwise  vorticity  is  defined  by 


_  Hi 

3y  3z 


Operating  on  w^  and  v^  results  in  the  desired  expression  for 


n.,(x,y,z,t)  -  l'  (x,z,t)  Tn(Y) 
N=*0 


(4.31) 


(4.32) 


(C) [Hq1 


«,„>  •  m>nu>,uD;  ic«2i  -  mi  iot*  1  -■=.  * 


o+8 


(4.33) 


where  D1>G1,  and  G2  are  defined  in  the  Appendix  A,  and  D2,  HQ,  C,  and  H2,  defined 
previously,  are  all  functions  of  a,B,  and  u.  This  expression  is  not,  however,  as 
complicated  as  it  looks.  The  residues  from  the  boundary  terms  are  negligible  in 
comparison  to  the  residues  from  A(a,8,u)  “0  so  that  these  terms  may  be  neglected. 
Furthermore,  since  [C][Hg]/A  is  the  solution  matrix  for  Vj(x,y,z,t),  the  above 
equation  involves  merely  a  pre-multiplication  of  the  v1  result. 

This  completes  the  solution  for  the  excited  first  orderfield  in  the  boundary  layer 
expressed  in  terms  of  physical  variables.  In  the  next  section  the  statistics  of  this 
result  are  derived. 


'j  ’,  .V  .  y _ _  .  ...... 


»  >  -  -  ^ 

V.  •”  *  'V.\'  •  .\v. 


;  '  \  -  ’  NV-\  *V  ■'  V  V. \>Y'  ■ 


Section  5.  Statistics  of  the  First  Order  Solution 

The  main  concern,  in  the  present  analysis  of  the  effect  of  free  stream  disturbances 
on  the  laminar  boundary  layer,  is  the  situation  where  the  free  stream  disturbances  are 
random.  The  vector  {q^},  the  Chebyshev  coefficients  of  the  free  stream  transverse 
velocity,  is  a  random  vector.  The  complex  outer  product  of  this  vector  results  in  a 
matrix  of  terms  which  include  the  known  statistical  properties  of  the  free  stream 
turbulence. 

The  solutions  for  the  first  order  velocity,  pressure,  and  vorticity  field,  derived  in 
the  previous  section,  were  shown  to  be  expressed  in  terms  of  the  vector  {q^} .  To  arrive 
at  the  statistical  properties  of  the  resulting  field,  the  complex  outer  product,  of  the 
Chebyshev  vectors  of  the  various  dependent  variables,  is  taken,  resulting  in  a 
relationship  between  the  statistical  response  in  the  boundary  layer  and  the  statistics  of 
the  free  stream  turbulence. 

In  what  follows  the  multidimensional  two-point  autocorrelation  functions  for  the 
three  velocity  components  are  derived.  The  derivation  of  the  statistics  of  the  transverse 
velocity  is  given  in  some  detail,  and  the  results  for  the  remainder  of  the  field  follows 
in  a  similar  manner,  and  therefore  the  results  only  are  given.  The  cross-correlations  for 
the  field  variables  can  also  be  derived  in  a  similar  manner,  and  this  is 
straightforward.  The  autocorrelation  for  the  streamwise  vorticity  is  also  given.  The 
other  vorticity  correlations  can  be  derived  in  a  similar  manner. 

The  solution  for  the  transverse  velocity,  v^,  has  the  Chebyshev  expansion, 

N 

v, (x,y,z,t)  =  I'  v1n(x,z,t)  Tn(Y)  (5.1) 

N"0 

The  two-point  autocorrelation  of  this  function  is  simply  <v^*,v^>.  In  terms  of  the 
Chebyshev  series  this  becomes 


<v1*,v1>  =  <v1*(x1  ,yi  ,z,t),v1(ic2,y2,z+5,t+T)> 

N  N 

=  l’  y*  <v1m*(x1,z,t),v1n(x2,z+i;,t+T)>Tn(Y1  )Tm(Y2)  (5.2) 

n=0  m=0 
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The  autocorrelation  function  is  therefore  a  double  Chebyshev  series,  the  coefficients  of 


2 

which  form  a  matrix  with  (N  +  1)  entries.  This  matrix  is  constructed  by  taking  the 

complex  outer  product  of  the  vector  {v1n}.  This  matrix  will  be  denoted  by  a  capital 
letter,  i.e.  [V,J.  Thus  [ V ^ ]  may  be  expressed  as 

(V1]  ’  {v1n*}  *  {v1n}T  (5*3) 

where  the  superscript  T  denotes  transpose.  The  vector  {v,n}  is  9iven  in  equation 
(4.16).  Before  proceeding  further,  it  is  important  to  notice  the  following.  The 
matrices  [I]  and  [G 1 ]  in  the  integrand  are  exponentially  smaller  than  the  main  term, 
as  will  be  shown  in  the  next  section,  and  may  therefore  be  neglected.  Then  substitution 
of  (4.16)  into  (5.3)  results  in 

tv,]  =  ////// 

-ffl 

(C*(a  ,6  ,ai  )  ]  [H*  ( a  ,6  ,ai  )]{q*Mq  }T[H  (a  ,8  ,u  >)T[C(a  ,B  ,a>)]  T 

_ 112  0111  n  m  0222 _ 2  2  2  _ 

A*(a,  ,8,  ,<0,  )4(a2,62,ui2) 

i[-a  x  -8  z+ui  t+a  x  +B,(z+?  )-u  (t+t)) 

e  *  da1da2dB1dB2dui1dm2  (5.4) 


However,  it  was  shown  in  Section  3  that 


“  [R]Q(-o1,a2,81,u1)S(»2-«1)B(B2-B1> 


(5.5) 


where  the  entries  of  the  matrix  (Rj  are  defined  in  equation  (3.17),  and  the  scalar  valued 
function  Q  is  defined  in  equation  (3.11).  The  Orr-Sommerf eld  operator  also  contains 
certain  useful  relationships  between  itself  and  its  complex  conjugate. 


D  *(a,B,w)  «  D^-a.B.-w) 

4  4 


D. (o,-B, w)  -  D  ( a, 6 ,ui) 
4  4 


a(8,-u)  ■  -a*(B,u) 


(5.6a) 

(5.6b) 

(5.6c) 


Proof  of  these  relations  can  be  found  in  the  appendix  of  the  paper  by  Tam  and  Chen 
(1976).  In  the  same  manner, it  can  be  shown  that 


H*( a, 8 ,<d)  -  HQ( -a, 8 ,ui) 


(5.6d) 


Substituting  (5.5)  and  (5.6)  into  (5.4)  results  in 

«  (C( a. ,8, -ai))  [H  (0,,8,-w))  [R]  [Hft(a,,6,(D)]T(C(a,,6,u)]T 

[V  )  -  ffff  - - 2 — 1 - °_i - 2 _ 

‘  1J  J1JJ  A(a1  ,8,-ui)4(a2,8,ai) 

»  ifa^^OjXj+Sc-uiT) 

Q(o1  ,a2>8,ui)e  do^do^dSdw  (5.7) 

This  integral  can  be  inverted  using  residue  theory.  Computational  aspects  of  this 
procedure  are  discussed  in  the  next  section. 

Since  the  procedure  is  similar  for  the  other  correlations,  only  the  necessary  details 
are  repeated  in  the  other  derivations. 


The  two-point  autocorrelation  for  u«  has  the  expression 
N  N 

<u* , u1 >  *  l'  l'  <u1n*(x1 ,z,t) ,u1m(x2,z+C,t+T)>T  (Y1 )Tm(Y2)  (5.8) 

n-0  m-0 

The  necessary  matrix  of  coefficients  is  given  by  the  complex  outer  product  of  the  vector 
(u1n)  given  in  equation  (4.28).  Neglecting  the  higher  order  terms  the  result  for  [Uj] 


[U1J  -  »,-«))  IV^o^Oj.S,!!!))  [H^a^Biw)  JT 


-1  T 

[o2  (a2,S,o>)]  •  Q(a1  ,a2 ,8,u’)e 


i  { a^x  1  +a2x2+8t-ii)T ) 


da^ajdSdu 


[V1  (a1  ,a2 ,8,<ii)) 


[HQ(a1 ,8,-u))  [R]  lHQ(  a2 , 8 ,u>)  ] T  [C(<*2,8  ,<■>)  )T 
i(a1 ,8,-(D)4(a2,8,w) 


(5.10) 


The  two-point  autocorrelation  for  w^  has  the  expression 


N  N 

<w1*.w1>  “  I'  I'  <%1*<Vz't)'^m(VZ+C't+T)>VVVV  <5-"> 

n«0  o-O 

with  the  necessary  matrix  of  coefficients,  [W^],  given  by  the  complex  outer  product  of 
the  vector  {w.^},  9lven  in  equation  (4.30).  Neglecting  the  higher  order  terms,  the 
resulting  expression  for  [W^]  is 

[w^  “  ////  82[D"1(a1,0,-w)nH2(a1,B,-<D)nvi(a1,a2,6,i.i)]  [H2(a2,B,u)]T 

—00  ** 

T  *  i(a1x1+a2x2+0c-uiT) 

[°2  *  Q<ai  >a2'8'“'e  da1da2dBdu  (5.12) 

where  [V^ )  is  given  in  (5.10). 

The  vorticity  correlations  can  also  be  derived.  For  example,  the  streamwise 
vorticity,  which  was  derived  in  the  previous  section,  has  the  following  two-point 
correlation, 

N  N 

<(!*,(),>  -  l>  l'  <ni*(x1,z,t),nim(x2,z+C,t+T)>Tn(Y1)Tm(Y2)  (5.13) 

n=*0  m=0 

Defining  the  associated  matrix  of  coefficients  by  [A^l,  after  neglecting  the  higher 
order  terms,  the  necessary  matrix  of  coefficients  is 

to  A 

(A,)  -  ////  B2(  (D1)  [D^cy  8, -ui)]  (H2(a1,B,-u))J-(I)]  (V1(a1,a2,B,u))  • 

—  00 

( (D^  (D21(a2,B,m)]  [H2(a2,8,u)]-[I]  1  TQ  (a1  ,a2,6 ,u)  • 

i  ( a  .x  1+a2x2-*-BC-o)T ) 

e  da^OjdSdu  (5.14) 

where  ( V ^ ]  is  defined  in  equation  (5.10). 

These  integrals  can  all  be  evaluated  in  a  similar  manner  using  residue  theory  in 
conjunction  with  numerical  integration.  The  calculation  centers  on  the  computation  of 


Section  6.  Inversion  of  the  Fourier  Integrals 

The  integrals  in  the  previous  section  are  of  the  type 


r(x,z,t)  =  Iff  ^4^4  ei( C,X+®Z~“t)dadBdw 
„  4(a,8,0)) 


where  F(a,8,io)  is  an  entire  function,  and  4(a,8,io)  is  a  determinant  expressing  the 
relationship  between  a  and  uj  for  particular  values  of  8,  often  referred  to  as  the 
dispersion  relationship.  Although  a  Fourier  transform  has  been  used  in  time  in  the 
analysis,  and  a  Fourier  transform  is  to  be  used  in  the  eventual  computation,  the  formalism 
for  inverting  the  transforms  is  derived  using  a  Laplace  transform  in  time.  The  formalism 
can  be  found  in  Chapter  2  of  the  monograph  by  Briggs  (1964).  Neglecting  8  for  brevity, 
the  expression  (6.1)  after  using  residue  theory  is 


M  «»tiCT  F(ot.  ,w) 

F(x,t)  =  2iri  J  f  — - - = - - 

k=1  '~*i0  ^  A(a'“>  a“a. 


i  ( a,  x-ut ) 
k 

e  du 


The  co-integration  is  deformed  down  to  the  real  axis  while  avoiding  all  real  values  of  a 
which  produce  complex  values  of  ui  with  positive  imaginary  parts.  In  addition  in  the  sum 
in  (6.2)  only  the  first  term  will  be  retained  as  it  produces  the  dominant  contribution.  A 
typical  spectrum  for  the  determinant  is  shown  in  Figures  (la)  and  (1b).  There  is  an 
isolated  unstable  mode  and  then  a  set  of  discrete  approximations  to  the  continuous 
spectrum  which  contribute  very  little  to  the  solution. 

Reinstating  8,  the  integral  (6.1)  can  now  be  evaluated  further.  First  the  ranges 
of  integration  for  w  and  8  may  be  simplified.  By  virtue  of  equations  (5.6),  the 
integrands  are  even  functions  of  8.  Using  the  [V.,]  matrix  as  an  example,  the  result  is 

00  CO 

(Vt )  *  2  /  Iff  I(  a1  ,a2 ,8,u)da1  da2du)d8  (6.3) 
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Figure  la  The  eigenvalue  spectrum  for  the  Blasius  boundary  layer  for  the 
case  of  spatial  stability  plotted  in  C-space  for  R  =  600,  u i  -  0.12,  0-0,  and 
.V  =  36. 
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Figure  lb  The  eigenvalue  spectrum  for  the  Blasius  boundary  layer  plotted 
in  o-spare  for  u.'  0.12.  3  0.,  R  -  600,  and  ,V  =  36. 
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i(a1 , 8 #  —to )  A <  c»2 ,  B ,  co ) 


Q(a1  ,a2,B,(u)e 


i (a1+a2>x 


(6.4) 


but  splitting  the  w  integral  into  an  integral  from  -■  to  0  plus  an  integral  from  0 
to  ao,  and  noting  that 

I(a1#a2#  8,-w)  *  I*( 8,a>)  (6.5) 

the  integral  becomes 
00  00 

[V1  ]  =2  ////  [I(a1,a2,B,u))+I*(-a1.-a2,B,(D)]da1da2dBda>  (6.6) 

0  —00 

The  result  (6.6)  is  proved  using  the  relations  given  in  equations  (5.6).  The  other 
integrals  can  be  treated  similarily. 

As  suggested  in  the  last  section,  the  Fourier  integrals  are  inverted  using  residue 
theory  in  and  a2<  while  the  cii  and  $  integrations  are  carried  out  using 

numerical  integration.  It  was  also  suggested  that  for  each  value  of  m  and  $»  there  is 
but  one  residue  in  each  of  the  a-planes  which  need  be  evaluated.  It  is  the  purpose  here 
to  see  what  other  poles  are  present,  illustrating  the  complete  solution  possibilities,  and 
where  and  why  many  of  these  poles  can  be  neglected.  The  form  for  the  solution  can  be 
expressed  generically  as 

<p(x)  =  /  K(  n,  X )  f(Sl)  dfl  (6.7) 

n 

where  <p  can  be  a  correlation  function  for  any  of  the  sought  field  variables,  and  X, 
representing  spatial  and  temporal  separation,  is  the  solution  space.  f(f!)  is  a  spectral 
function  representative  of  the  random  forcing  function,  defined  in  a  spectral  volume 
fi.  K(fi,X)  is  a  deterministic  transfer  function  which  represents  the  operator.  This  is 
the  classical  form  for  a  random  forced  vibration  problem.  There  are  obviously  two  sources 
for  poles  in  the  integrand,  the  transfer  function  K,  and  the  forcing  function  f. 
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To  illustrate  these  ideas  it  is  necessary  to  examine  only  the  expression  for  the 
Vj  autocorrelation  and  the  u^  autocorrelation.  The  others  are  derived  from  these  or 
have  the  same  poles. 

Equation  (5.7)  of  the  last  section  defines  the  matrix  necessary  for  the  v^ 
correlation.  The  matrix  (R)  and  the  function  Q  are  representative  of  the  forcing 
function.  Because  the  expression  is  in  terms  of  matrices,  the  matrix  [R]  is  sandwiched 
in  the  midst  of  the  transfer  function.  The  transfer  function  is  represented  by  the 
remaining  matrices,  [C]  and  (HQ] ,  and  the  determinant,  which  appears  in  the 
denominator. 

The  forcing  function  has  one  pole.  If  the  integration  over  is  performed  first, 

there  is  a  pole  at  a  »  iy  -  o^,  as  can  be  seen  in  equation  (3.11).  However,  when  this 
pole  is  substituted  into  the  exponential  terms,  explia^x^,  the  result  is 
exp(-yx2  -  ia.jX2).  Since  y  >  0,  and  ^  ia  proportional  to  the  Reynolds  number,  this 
term  is  exponentially  small  in  comparison  to  contributory  poles  with  positive 
exponentials.  Therefore,  this  pole  can  be  safely  excluded  from  the  calculation.  If  this 
pole  had  a  more  significant  contribution  it  would  result  in  the  forcing  function  imposing 
its  character  on  the  boundary  layer.  This  is  a  phenomenon  which  is  often  noticed  in 
experiments  with  large  amplitude  free  stream  tubulence.  Reshotko  (1976)  refers  to  this 
phenomenon  as  "high-intensity  bypass”.  It  is  possible  that,  for  the  case  of  strong  free 
stream  turbulence,  this  pole  dominates  the  other  poles. 

In  the  solution  for  the  v^  autocorrelation  matrix,  the  remaining  contributing  poles 
are  due  to  the  roots  of  A(a,B,u>)  m  0.  These  are  the  natural  modes  of  the  Orr-Sommerf eld 
operator,  discussed  in  some  detail  by  Bridges  and  Morris  (1985).  There  will  be  both 
upstream  traveling  modes  (Re(a)  <  0),  and  downstream  traveling  modes  (Re(a)  >  0). 
Numerical  evidence  suggests  that  the  modes  with  Re (a)  <  0  form  a  continuous  spectrum  and 
are  all  damped.  Whether  there  are  discrete  upstream  modes  or  unstable  modes  is  unknown. 
The  remaining  spectrum  can  be  broken  up,  for  discussion  purposes,  into  three  parts;  (a) 
the  single  dominant  mode,  (b)  the  higher  modes,  and  (c)  the  continuous  spectrum  (or 
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discrete  approximation  of  it).  All  three  contribute  poles  to  the  solution.  However,  in 
the  range  of  integration,  there  is  a  band  of  frequencies  and  spanwise  wavenumbers  which 
produce  eigenvalues  with  negative  imaginary  parts.  The  other  poles  always  have  negative 
exponentials.  Therefore,  the  single  dominant  mode  produces  a  significantly  greater 
contribution,  and  the  other  poles  can  be  safely  neglected.  A  similar  discussion  can  be 
put  forward  for  the  auto-correlation.  The  matrix  necessry  for  this  correlation  is 

given  in  equation  (5.9)  and  (5.10).  Since  the  v1  autocorrelation  matrix  appears  in  this 
expression,  it  will  have  the  poles  discussed  above,  with  the  same  argument  holding  for 
retention  of  the  dominant  mode  only.  However,  in  this  expression,  there  are  two 
additional  sources  for  poles.  The  operator  matrix  Djta.S.iu)  has  poles.  These,  however, 
can  be  safely  neglected.  In  Appendix  C  it  is  proved  that  the  poles  of  this  operator  are 
heavily  damped  and  offer  little  to  the  solution.  The  other  source  is  the  matrix 

H^a.B.ui).  This  matrix  is  a  discretization  of  the  right  hand  side  of  equation  (4.23). 

2  2 

Note  that  there  is  a  term  a  +  B  ,  in  the  denominator.  This  gives  rise  to  poles  at 
a  =  +iB.  Since  there  is  no  branch  cut  in  the  a-plane  between  a  ■  -iB  and  a  *  +iB, 
the  contour  of  integration  is  deformed  such  that  only  the  pole  at  a  *  +ifi  will  be  inside 
of  the  contour.  But  this  pole  produces  an  exponential  weighting  of  the  form  exp(-Bx), 
again  an  exponentially  small  term. 

In  conclusion,  it  is  clear  that  the  forcing  function  and  operator,  or  transfer 
function,  are  rich  in  their  possibilities  for  solution,  but  the  exponential  term  in  the 
integrand  filters  out  many  of  these  possibilities  leaving  the  dominant  mode  as  the 
residuum. 

With  these  ideas  in  mind  the  Fourier  transforms  can  now  be  inverted.  Using  the 
[v^  autocorrelation  matrix  for  illustration  and  setting  £  •  t  =  0  for  brevity,  the 
result  after  integration  over  and  a2  is 

ao 

[V  ]  =  -8t2  //  (1(8, co)  +  I*  (  8  ,w) )  dudB  (6.8) 

1  0 
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where 


1(8, w) 


[C(-a*,8»-u)]  [H0(-a*,8,-w)]  iR]  [H0(ak,8,uD]T[C(ak,8,u)))T 


A(,>(-a*,e,-w)A<1)(V8'“) 

1 

Qt-ajJ.a^S.wJe 


i(a  -of) 
k  k 


(6.9) 


and  ak  satisfies  Afo^jS.aj)  m  0,  and 

„  A0*xlz*t  f  1  i 

Q(-o*,ak,6.U) - —  [~  7  t(Vaj>)]« 


2 

-8  A 


2* 


where 


02  ■  *xV+  ‘z2®2  +  <“  -  VV 


and 


(6. 10a) 


(6.10b) 


A(1)  =  dA/da  . 

Before  proceeding  further,  it  is  obvious  that  it  will  be  necessary  to  evaluate  the 
derivative  of  a  determinant  at  a  singular  point;  A(ak,8,u)  “  0,  where  it  is  assumed  that 
ok  is  a  simple  root. 

A  general  efficient  method  for  evaluating  the  local  derivative  of  a  determinant  at 
singular  and  non-singular  points,  as  well  as  higher  derivatives,  has  recently  been 
developed  by  Bridges  and  Vaserstein  (1985).  A  brief  outline  of  the  necessary  details  are 
presented  in  Appendix  B. 

After  the  residue  theory  is  used  in  the  integral  over  o1  and  o2,  the  remaining 
work  is  a  numerical  integration  of  an  integral  of  the  form 


8  <D 

F  «  /  /  f(6,w)dwd8  (6.11) 

0  0 

and  this  is  carried  out,  in  the  classical  manner,  using  Simpson's  rule. 

This  completes  the  computational  details  needed  to  evaluate  explicitly  the 
statistical  properties  of  the  excited  solution  in  the  boundary  layer.  In  the  next  section 
results  of  the  evaluations  are  presented. 
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Section  7.  Results 


Since  the  zeroth  order  solution  is  the  well  known  Blasius  solution,  the  present  study 
need  only  concentrate  on  the  perturbation  to  this)  the  statistical  results  for  the  first 
order  field.  The  set  of  parameters  governing  the  first  order  set  of  solutions  is 
numerous.  Representative  of  the  free  stream  turbulence  are  the  length  scales 
tx,  ly,  and  t'.e  time  scale,  t 1,  and  the  decay  rate  y.  The  statistical  response 

can  be  evaluated  in  a  variety  of  ways.  The  response  consists  of  auto-  and  cross¬ 
correlations  (or  the  spectra)  of  the  three  velocity  components,  the  pressure,  and  the 
three  components  of  vorticity.  As  the  correlations  and  spectra  are  interchangeable,  only 
the  correlation  results  are  considered.  Each  of  the  correlation  functions  are  dependent 
upon  the  separation  distances  in  x,  y,  z,  and  t. 

There  is  also  numerical  accuracy  to  consider.  The  numerical  considerations  are  (a) 
the  order  of  the  series  expansions  in  terms  of  Chebyshev  polynomials,  and  (b)  the 
numerical  integration  in  m  -  8  space.  In  all  the  results  presented,  N,  the  order  of 
the  Chebyshev  expansion  is  greater  than  or  equal  to  24.  Numerical  tests  with  other  values 
of  N  have  shown  this  to  be  sufficient.  For  the  numerical  integration  of  the  spectral 
integrands  in  w  -  8  space,  after  testing  several  nets,  a  value  of  11  points  in  u  and  9 
points  in  8  was  chosen. 

Rather  than  concentrate  on  a  parametric  study,  the  main  theme  will  be  to  elucidate 
the  nature  of  the  three  dimensionality  in  the  response  of  the  boundary  layer.  For  the  set 
of  results  presented,  the  following  parmeters  are  used  for  the  free  stream  turbulence; 

=  1/5,  =  1/5,  =  1/5,  =  1/5,  and  y  =  .01,  where  these  parameters  have  been 

normalized  by  the  Blasius  boundary  layer  thickness  and  the  free  stream  velocity. 

The  subscripts  for  the  first  order  solutions  will  be  dropped  for  brevity.  The  first 
order  y-velocity  will  be  v,  etc. 

Figures  2,  3,  and  4  show  the  transverse  distributions  of  the  excited  response  at  R  = 
1000  (R  here  is  based  on  displacement  thickness).  In  Figure  2,  the  rms  distribution  for 


the  three  velocity  components  is  plotted  (this  is  the  square  root  of  the  autocorrelation 
function  at  zero  delay)  versus  the  distance  across  the  boundary  layer  (y/6  =1  is  at  the 
edge  of  the  Blasius  boundary  layer).  The  most  important  feature  here  is  that  all  three 
components  are  of  a  similar  magnitude,  emphasizing  the  three-dimensionality  of  the 
response.  The  magnitudes  have  been  scaled  such  that  the  maximum  of  the  streamwise 
component  is  unity  with  the  other  components  scaled  accordingly.  The  actual  magnitude 
vary  widely  as  the  parameters  governing  the  free  stream  turbulence  are  varied  (the  excited 
magnitudes  range  from  order  10  to  order  1000).  However,  the  qualitative  features  of  the 
excited  field  are  invariant  to  the  changes  in  magnitude.  Therefore  in  these  results  the 
emphasis  will  be  on  the  qualitative  features  and  comment  on  the  relationship  between  the 
free  stream  parameters  and  the  excited  magnitudes  must  be  reserved  until  a  more  complete 
parametric  study  can  be  performed. 

It  is  also  interesting  to  note  the  hump  in  the  distribution  of  urms  near  the 
wall.  This  is  a  distinctive  feature  of  the  exerimental  results  of  Klebanoff,  et.  al 
(1962)  and  will  be  discussed  further  at  the  end  of  this  section.  The  maximum  of  each  of 
the  functions  occurs  near  the  wall  with  the  exception  of  vms-  This  in  fact  holds  true 
for  other  Reynolds  numbers  (up  to  3000)  as  well.  The  maximum  of  Vrm8  occurs  at  about 
60%  of  the  boundary  layer  thickness.  According  to  Klebanoff,  et.  al.  (1962)  this  is  the 
region  where  non-linear  breakdown  occurs.  It  is  likely  that  an  instantaneous  shear  layer 
is  set  up  in  this  region  due  to  the  fact  that  3v/3y,  on  average,  is  zero  there,  which  is 
highly  unstable. 

In  Figure  2  the  relative  magnitudes  are  <vrms)max  =  .38  (urms)max  and  ( *»rms ) = 

.77(unns)(nax.  In  Figure  3  the  rms  value  of  the  streamwise  and  spanwise  vorticity  are 

plotted  across  the  boundary  layer.  Here  ( Hx  )  =  .77 ( Dz  )  .  The  streamwise 

‘  rms  max  rms  max 

vorticity  is  very  close  in  magnitude  to  the  spanwise  vorticity  emphasizing  the  three 
dimensional  nature.  This  also  shows  that  streamwise  vorticity  can  be  excited  by  the 
presence  of  infinitesimal  free-stream  turbulence  whereas  it  was  previously  thought  that  it 
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was  due  to  a  secondary  internal  instability.  In  Figure  4  the  distribution  (in  the 
computational  domain)  of  the  Reynolds  stress  is  plotted.  The  left  axis  is  the  computation 
domain  and  the  right  axis  is  the  Blasius  variable  (y  =  5  is  the  edge  of  the  Blasius 
boundary  layer).  The  maximum  occurs  at  about  20%  of  the  boundary  layer  thickness.  The 
qualitative  features  shown  here  vary  little  as  the  Reynolds  number  is  increased  to  2000 
and  3000.  The  function  governing  the  Reynolds  stress  gets  a  little  broader  and  the  ratio 
of  the  streamwise  to  spanwise  vorticity  maximums  drops  to  .58  at  R  =  3000  (compared  with 
.77  at  R  *  1000).  Otherwise  the  results  are  similar. 

However  these  results  are  for  zero  separation  in  space  and  time.  In  keeping  with  the 
theme  of  elucidating  the  three-dimensional  effects,  results  for  non-zero  spanwise 
separation  will  be  shown.  In  Figure  5  the  normalized  <u,u>  velocity  auto-correlation  is 
shown  as  a  function  of  spanwise  separation.  This  is  for  a  fixed  point  in  the  boundary 
layer  (a  point  in  y  which  is  at  26%  of  the  boundary  layer  thickness).  The  results  for 
R  «  1000,  2000,  and  3000  are  included.  5  is  normalized  by  the  Blasius  variable, 
therefore  C  ”  5  would  correspond  to  a  boundary  layer  thickness.  The  result  shows  a 
definite  finite  structure  in  the  spanwise  direction  with  increasing  spanwise  length  scale 
with  increasing  Reynolds  number.  This  characteristic  scale  is  not  a  function  of  the  much 
smaller  length  scale  of  the  freestream  turbulence  but  is  governed  by  the  unstable  response 
of  the  boundary  layer  over  a  wide,  but  band-limited,  range  of  three-dimensional  normal 
modes . 

In  Figure  6  a  similar  result  for  the  streamwise  vorticity  is  plotted  at  the  same 
location  in  y  as  the  <u,u>  correlation.  Again  there  is  evidence  of  a  three-dimensional 
character  indicative  of  streamwise  vortices  of  large  scale  excited  in  the  boundary 
layer.  The  typical  length  scale  in  the  spanwise  direction  of  the  streamwise  vorticity  is 
also  growing  with  increasing  Reynolds  number. 

The  main  new  feature  of  the  results  shown  so  far  is  the  prevalent  three 
dimensionality.  Since  the  formulation  of  the  problem  is  for  free  stream  turbulence  of 
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infinitesimal  amplitude,  this  suggests  that  in  any  real  situation  a  three  dimensional 
response  of  the  nature  shown  will  be  excited.  However  it  is  instructive  to  ask  how  the 
solution  might  differ  if  the  free  stream  disturbances  were  strictly  two-dimensional.  In 
this  case  the  spectral  integrands  contain  a  Dirac  delta  function  in  B.  Therefore  only 
the  8  =  0  point  of  the  integrand  is  taken  resulting  in  a  two-dimensional  response  in  the 
boundary  layer.  The  evaluation  of  the  normalized  <v,v>  and  <u,u>  auto-correlations, 
for  zero  delay,  at  R  =  1000,  for  the  case  of  a  two-dimensional  free  stream  are  shown  in 
Figure  7  where  the  distribution  of  the  rms  u  velocity  is  plotted.  This  u  result  is 
quite  different  from  the  three-dimensional  result.  The  2-D  result  is  much  smoother  and 
doesn't  contain  the  extra  "hump"  near  the  wall  which  is  prevalent  in  the  3-D  results. 

This  result  can  be  compared  with  the  experimental  measurements  of  Klebanoff,  et  al. 

(1962).  In  their  experiments  on  natural  transition  they  measured  the  rms  distribution  of 
the  u  velocity.  Figure  8  is  a  reproduction  from  their  paper.  This  can  be  compared  with 
the  three  dimensional  results  in  Figure  2  and  the  2-D  result  in  Figure  7.  Although  the 
result  of  Klebanoff,  et  al.  is  for  a  Reynolds  number  of  2800,  Figures  2  suggests  that  the 
structure  of  the  rms  distribution  varies  little  with  Reynolds  number.  The  presence  of  the 
extra  "hump"  in  the  result  of  Klebanoff,  et  al.  suggests  that  there  is  a  3-D  excitation 
present  in  their  experiment.  This  structure  is  not  unique  to  this  particular  figure  of 
Klebanoff,  et  al.  either.  Most  of  the  experimental  results  displayed  in  the  paper  of 
Klebanof,  et  al.  have  a  similar  structure. 
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Section  8.  Discussion 

Using  the  classical  theory  of  random  forcing  of  a  linear  MDOF  vibrating  system,  the 
random  excitation  of  the  Blasius  boundary  layer  has  been  analyzed.  The  Blasius  boundary 
layer  is  treated  as  a  MDOF  vibrating  system  subjected  to  a  random  forcing  functions  the 
free  stream  turbulence.  The  main  result  of  the  present  analysis  is  the  prevalence  of 
three  dimensionality  in  the  result.  This  suggests  that  boundary  layers  will  have  three- 
dimensional  velocity  and  vorticity  fields  excited  in  them  due  to  the  presence  of 
infinitesimal  free  stream  turbulence  (a  virtual  surety) .  Furthermore  the  three 
dimensionality  is  a  linear  phenomena,  and  can  therefore  be  referred  to  as  a  primary 
instability  rather  than  a  secondary  instability. 

There  does,  however,  remain  some  unclear  points.  The  most  important  of  these  is  what 

can  be  referred  to  as  the  paradox  of  boundary  layer  transition.  This  can  be  explained  as 

follows.  In  Section  2  it  was  pointed  out  that  for  infinitesimal  amplitude  free  stream 

2 

turbulence,  the  terms  of  order  e  can  be  dropped,  and  then  the  governing  equations  for 
the  perturbations  arising  in  the  boundary  layer  are  linear.  Talcing  the  limit  of  these 
equations  as  y  ♦  «  should  result  in  the  governing  equations  for  the  free  stream 
disturbances.  Unfortunately  the  eigenvalues  of  this  set  of  equations  form  a  continuous 
spectrum  and  the  eigenfunctions  are  harmonic.  Since  any  set  of  disturbances  will 
essentially  be  expanded  in  terms  of  these  eigenfunctions,  this  set  of  equations  will  not 
support  disturbances  with  discrete  eigenvalues.  Therefore  this  free  stream  turbulence 
will  not  be  capable  of  exciting  Tollmien-Schlichting  waves  in  the  boundary  layer. 

However,  it  is  generally  accepted  that  Tollmien-Schlichting  waves  do  actually  exist  in 
experiments,  and  that  linear  theories  are  in  agreement  with  this.  This  is  the  paradox  of 
boundary  layer  transition. 

This  was  overcome  in  the  present  study  by  assuming  that  the  flat  plate  was  situated 
in  a  large  wind  tunnel  and  therefore  the  set  of  equations  in  the  free  stream  is 
approximately  correct.  However,  in  a  more  complete  analysis,  it  may  be  necessary  to 
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A  similar  result  occurs  when  f(t)  is  a  random  function.  If  f ( t )  is  random  and 
stationary  in  time,  a  solution  cannot  be  found  because  the  transfer  function  has  a  double 
pole  on  the  real  uj  axis.  However,  if  f(t)  is  random  and  non-stationary  (decaying  in 
t,  for  example)  then  the  integrals  exist  and  a  statistical  solution  may  be  found. 

This  is  essentially  what  happens  in  the  boundary  layer  analysis.  The  neutral  curve 
is  a  singularity  of  the  linear  theory  if  the  free  stream  disturbances  are  homogeneous  in 
the  streamwise  direction.  However  the  governing  equations  for  the  free  stream  turbulence 
have  a  parabolic  term,  proportional  to  1/R,  which  suggests  a  variation  in  amplitude  of  the 
solution,  and  additionally,  on  physical  grounds  it  is  expected  that  free  stream 
turbulence,  in  the  absence  of  a  mean  shear,  will  decay.  Therefore  the  use  of  a  set  of 
disturbances  with  decaying  amplitude  is  the  proper  choice.  But  this  nevertheless  leads  to 
speculation  on  the  effect  of  non-linearity  on  the  singularity. 

The  SDOF  equation  can  again  be  used  as  an  analogy. 


+  <i)„ 


ef(t) 


(8.3) 


This  is  the  well  known  undamped  Duffing  equation  (Nayfeh,  Chapter  9).  To  find  solutions 

2/3 

for  w  near  (dq  one  defines  a  detuning  parameter  u>  “  a)Q  +  e  a.  Then  it  can  be 
formally  shown  that  <p  must  be  expanded  in  a  series 


?(t)  »  l  t  <p  (t) 

*  n 
n=1 


(8.4) 


This  shows  that  a  forcing  function  of  order  e  excites  a  response  of  order 
e  ^  (as  €  *  0)  near  resonance.  A  similar  argument  can  be  used  for  the  boundary 
layer.  A  standard  Stuart-Watson  expansion  for  the  non-linear  modes  near  the  neutral  curve 
suggests  that  the  amplitude  correction  to  the  linear  eigenvalues  is  of  the  order  of  the 
amplitude  squared.  Therefore  in  order  to  bring  in  the  non-linear  terms  correctly  to 
overcome  the  singularity,  the  solution  in  the  boundary  layer  should  be  expanded  in  powers 
of  E  ^  where  e  is  the  amplitude  of  the  free  stream  disturbances.  The  main 
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implication  ia  that  disturbances  in  the  free  atream  of  order  e,  which  are  hoetogeneous 
in  x,  with  diaperaion  propertiea  near  thoae  of  the  neutral  curve  will  excite  a  reaponae 
in  the  boundary  layer  of  order  e1^  (as  e  ♦  0).  The  details  of  such  a  non-linear 
analysis  are  not  given  here  as  the  complete  analysis  will  be  the  subject  of  a  future 
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Appendix  A 

The  matrix  [G|]t  first  used  in  equation  (4.14)  can  be  derived  in  the  following 
manner.  It  is  a  result  of  the  Chebyshev  discretization  of  the  boundary  condition 
function  g.,  (x,y,z,t) . 

If 

y  -  2(1  +  Y)/(1  -  Y)  (A. 1 ) 


and 


and 


then 


exp(-y)  »  l'  d  T  (Y) 


n“0 


yexp(-y)  *=  J'  e  T  (Y) 


n»0 


N 


g1n(x,z,t)  -  l'  (-1)3q3  (x,z,t)  (d  +  e  (1  -  j^)) 
j-0  n  n 


’In' 


therefore  the  matrix  [G.)  has  entries 


[G1]i,j  “  (_1,j  (di  +  ei(1  “  32)>  for  1  =  °<  N 

j  -  0,  1,  ...»  N 


(A. 2) 


(A. 3) 


(A. 4) 


(A. 5) 


The  matrix  [Gj]  arises  from  the  Chebyshev  discretization  of  the  boundary  condition 
function  g2(x,y,z,t)  and  it  has  entries 


[G2]ij  -  <-1)3  dt  for  i  -  0,  1,  ....  N 


and  j  ■  0,  1,  ,  N 


(A. 6) 


where  d^  is  given  in  (A.2). 

The  matrix  (D^)  is  a  discretization  of  the  operation  md/dY.  The  metric  m(Y)  has 
the  Chebyshev  expansion. 


m(Y)  -  j  mQ  +  m^T1  (Y)  m^iY) 


(A. 7) 


where  mQ  »  3/4,  m^  ~  -1/2,  and  m2  *  1/8.  If  we  define  a  pentadiagonal  matrix  [M] 
such  that 
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then  the  matrix  [D«]  has  entries 


tD1]i  -i  “  ^  M,  J1  -  (-D3+lt)j 

1,3  K=0  1,K 


for  i  •  0,1,  .  . . ,  N 


and  j  »  0,1,  . . . ,  N 


Appendix  B 


Consider  the  matrix  [D  (a)I,  which  depends  on  the  scalar  a  in  the  following  way. 


f  C.aP_j 


where  the  are  square  complex  matrices.  The  determinant  of  D^ta)  is  defined  as 


A(a),  and  its  derivative  is 


A(1)(a)  -  -r^  A( a) 
da 


It  is  necessary  in  the  present  work  to  evaluate  the  derivative  at  (simple)  latent  roots. 


where  A  (a.  )  «  0. 


The  classical  expression  for  the  derivative  of  a  determinant  can  be  constructed  in 


the  following  way.  Define  D.^(a)  as  the  jth  column  of  D^(a),  and  (a),  as  the 


derivative  of  Dp(a)  with  respect  to  a.  Then  if  N  is  the  order  of  the  matrix  D^(a), 


the  derivative  of  the  determinant  of  D  (a)  can  be  expressed  as 

P 


(  1 ) 

A  ( a)  =  V  A .  ( a) 


.  v-S-;-.-:.': 


where  A.,  (a)  is  the  determinant  of  a  matrix  whose  jth  column  is  D,j  (a)  and  the 
remaining  columns  are  those  of  D^fa).  Unfortunately  this  expression  requires  an 
operation  count  of  0(N4)  for  each  point  a^,  which  is  clearly  prohibitive. 

with  a  slight  rearrangement  of  the  matrices,  the  local  derivative  of  a  determinant 
can  be  evaluated  in  a  straightforward  and  efficient  manner.  The  idea  is  to  arrange 


Dp (a)  in  the  following  Taylor  form 


D  (a)  »  y  F  (a  -  a. 
P  1  k 

j»0 


1  dJ 

F,  *  — T"  — 7-  D  (a) 

!  .  i  .  j  p 

3  da  a-a 


However,  only  the  first  two  matrices  need  to  be  retained  to  determine  A  '  ( a  )  exactly. 

D  (a)  »  Fn  +  F, ( a  -  a.  >  +  ...  (B.6) 

p  U  1  K 

Obviously,  fq  is  singular  since  Alo^)  *  0,  therefore  FQ  can  be  brought  to  the  form 


A  “  Fo  “ 


0000---0 
Ox  x  x  x 
0  x  X  x 


The  first  row  and  first  column  of  Fq  are  brought  to  zero  and  the  remaining  N  -  1  by  N 
-  1  minor  is  brought  to  upper  triangular  form,  with  the  same  operations  performed  on 
F  ^ .  if  B»F1  (F1  after  operations),  then 

A(1,(ak)  =  (b1t)  •  (a22)(a33,  ...  (a^)  (B.8) 

The  derivative  of  the  determinant  is  the  product  of  the  determinant  of  the  principal  minor 
of  A  times  the  upper  left  entry  of  B. 

It  should  also  be  apparent  that  in  the  limit  as  a  ♦  o^,  there  is  only  one  non-zero 
entry  in  the  cofactor  matrix. 


limit  C(a) 
a  ♦  ok 
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where  =  (a22^a33^  *  •  •  *aNN^  i®  the  determinant  of  the  principal  minor  of  h.  In 

the  computation,  however,  (B.8)  appears  as  a  divisor  to  (B.9),  therefore  the  determinant 
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Theorem:  The  secondary  operator  L3  governs  the  streamwise  and  spanwise  velocity  pertur¬ 
bations  ir.  -he  Orr-Sommerfeld  problem. 


L2T 


d2  T 


dy2 


tR(uj  al’(y))  (oJ-dJ)]T 


(CM) 


where  u,3.  R  6  R  with  R  >  0,  a  G  C  is  the  “spatial”  eigenvalue 

If  T  is  square  integrable  in  y  f  [0,oo),  T (0)  =  0,  T(y)  0  sis  y  —  oc,  and  t/( y)  is  non¬ 
negative  and  bounded  for  all  y,  then  the  spatial  eigenvalues  of  La  are  stable,  i.e.  Im(a)  >  0, 
for  any  R  >  0 


Proof:  Set  LjT  -  0  and  multiply  thru  by  T",  the  complex  conjugate  of  T  and  integrate  over 

y 


/"“idT  2  r06  /*> 

1^;  dy*(i*;R-  a*  -  0*)  j  r^dy-iaRj  l'(y)|T  2dy  0  (C  2) 


With  C/(y)  non-negative  and  bounded  for  all  y  the  three  integrals  may  be  represented  by 
positive  definite  constants 


-a3  ■+•  (iujR  -  a2  -  02)a)  -  iaRa2  =  0 


(C.3) 


which  may  be  solved  fGr  a, 


iRai 

a  = - -  < 

2a, 


-°-{l  ±  t/ 1  +  ( 

1 1  ^  V  <12  R 

But  the  square  root  of  a  complex  number  is 


2a,0  2  40,03  ,4-ja^ 


a?R2  alR 


(CM) 


sjz  =  sjx  +  «y  =  ±{y-(|*l  +  x)  +  «'!si?n(y)'y/i(:*j  -  x)  } 


(C5) 


so  that 


2a  1 3  4aiv’3  t'4wa 

V,  +  (^r)  +af^' 


a2/? 


:---  a4  -  ra5 


(C  6) 


where  a4  and  as  are  positive  definite  constants  and  the  sign  of  u  is  ignored  as  it  has  no  effect 
on  the  outcome.  Then  from  (C  4) 


1  Ra2 


2a ) 


(1  ±  (a«  -  IQs)) 


(C 


Choosing  only  the  downstream  directed  wave  (Re(a)  >  0), 


Q‘=  2a7^,K  0  2a, 


(C  8) 
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Thus  the  eigenvalue  is  stable  is  a«  >  1  and  unstable  is  a(  <  1.  Now 


I  1,  1  ,2ai^u  4a,a3  1 

“«=lV*  2  >  "  n 2 W  2  a\W  \ 


ft.  /2ai^\2  40,03  l!  16w!aji! 

H,4(a2«)  +  aWi  +  aT«^/ 


Clearly  in  the  limit  as  R  —  oo,  a4  — *  1,  and  for  all  finite  values  of  R,  at  >  1.  This  proves  the 
theorem. 
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